Dark matter transport properties and rapidly rotating neutron stars 
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Neutron stars are attractive places to look for dark matter because their high densities allow 
repeated interactions. Weakly interacting massive particles (WIMPs) may scatter efficiently in the 
core or in the crust of a neutron star. In this paper we focus on WIMP contributions to transport 
properties, such as shear viscosity or thermal conductivity, because these can be greatly enhanced 
by long mean free paths. We speculate that WIMPs increase the shear viscosity of neutron star 
matter and help stabilize r-mode oscillations. These are collective oscillations where the restoring 
force is the Coriolis force. At present r-modes are thought to be unstable in many observed rapidly 
rotating stars. If WIMPs stabilize the r-modes, this would allow neutron stars to spin rapidly. This 
likely requires WIMP-nucleon cross sections near present experimental limits and an appropriate 
density of WIMPs in neutron stars. 

PACS numbers: 95.35.+d, 97.60.Jd, 26.60.Dd, 66.20.-d 
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Dark matter has proven to be elusive. Weakly interact- 
ing massive particles (WIMPs) have likely not yet been 
seen in the laboratory. In astrophysics, neutron stars 
may be promising places to look for WIMPs because 
they are the only known objects dense enough for WIMPs 
to interact repeatedly. WIMPs could condense to form 
a black hole and destroy their host neutron star [H-Q- 
Perez-Garcia et al. suggested that WIMPs could convert 
a (hadronic) neutron star into a strange star [6]. Al- 
ternatively, WIMP annihilation could heat neutron stars 
0, 0, HI- Finally Leung et al. considered the impact of 
dark matter on the frequencies of neutron star oscillation 
modes dGl]. 

In this paper we focus on WIMP contributions to 
transport properties. A few WIMPs may greatly increase 
the shear viscosity or thermal conductivity of neutron 
star matter because WIMPs have very long mean free 
paths Xw that efficiently transport momentum and en- 
ergy. For example Xw may be 10 10 times longer than 
the mean free path of an electron. This means that, per 
particle, a WIMP will contribute dramatically more to 
the shear viscosity. Note that by focussing on transport 
properties we are taking advantage of the long WIMP 
mean free path. This is the very thing that often makes 
studying dark matter so difficult. 

We start with a model where is smaller than the 
size of a neutron star. Next we calculate the WIMP 
contribution to the shear viscosity of dense neutron rich 
matter and speculate that this shear viscosity may damp 
r-modc oscillations in rapidly rotating neutron stars. Fi- 
nally we discuss WIMP contributions to the thermal con- 
ductivity. 

WIMPs could be trapped in the core or in the crust 
of a neutron star. We consider both possibilities in turn. 
First, let us assume Xw is dominated by scattering from 



nucleons in the core of a neutron star 
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with o Wn an effective transport cross section for WIMP- 
nucleon scattering that includes the effects of Pauli block- 
ing, see below. The nucleon density is n n . We assume 
n n 3> nw where nw is the WIMP density and we neglect 
WIMP self-interactions. We require that A^ is not much 
larger than the size of the star. Consider 
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10 km. 
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For simplicity we model a neutron star as a uniform den- 
sity sphere with say n n 2uq = 0.32 fm -3 . Here no 
is normal nuclear density. Note, the central density of a 
neutron star is expected to be significantly larger than 
no. Equation [2] requires 
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3 x 10 



" 45 cm 2 . 
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For Mw near 50 GeV and assuming spin independent 
interactions, this is slightly below - but near - the recent 
experimental limits of the CDMS and EDELWEISS [U 
and XENON 100 [l2j experiments. For larger or smaller 
masses, Eq. [3] is significantly below the experimental 
limits. 

Nucleon Pauli blocking may reduce the effective cross 
section awn in the star so that Xw is too large. In 
this case spin dependent interactions could still pro- 
duce the required mean free path, since they are much 
more weakly constrained by experiment. Alternatively, 
WIMPS could be trapped in the inner crust of the neu- 
tron star because of coherent WIMP-nucleus scattering. 
Although the density is lower in the crust, the coherent 
cross sections may be large. 

We now consider WIMP scattering in the neutron star 
inner crust. Coherent WIMP nucleus scattering cross 
sections are proportional to A 2 where A is the mass num- 
ber. We choose to normalize this cross section per nu- 
cleon. Therefore coherence increases the cross section per 
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nucleoli by a factor of A. This coherence is reduced by 
the form factor of the nucleus F(q) with q the momen- 
tum transfer and F(q = 0) = 1, F(q — > oo) = 0. In 
addition, in the dense medium other nuclei can screen 
the scattering. This is described by the static structure 
factor Sa (q)- The static structure factor adds coherently 
scattering from all of the different nuclei [13(. Screening 
ensures that Sa(q = 0) is small and the static structure 
factor is normalized Sa(q — » oo) = 1. 

The differential WIMP-nucleus scattering cross section 
(per nucleon) in the medium da/dil is 

da/dn^da/dn\oAS A (q)F(q). (4) 

Here da/dfl\o is the free space WIMP-nucleon scattering 
cross section. Integrating over scattering angles gives the 
total cross section 

a = J dndcr/dn = (Jo{A) (5) 

Here <r = / dflda/dfl\o and Eq. [5] defines the effective 
number of nucleons 

(A)=A(S A (q)F(q)) (6) 

In general (A) < A because SA(q) < 1 for small q and 
F(q) < 1. We expect the product SA(q)F(q) to peak 
at intermediate momentum transfers near q* « 0.3 fm -1 
[25j . If the momentum transfer range for the integral in 
Eq. [S]includes q* then (A) may only be somewhat smaller 
than A say ssl/2 |25j . However if the WIMP momentum 
is too small, so that q* is not part of the integration range, 
than (A) can be <C A because Sa{q) "C 1 for small q. 

We now estimate the optical depth r for WIMP scat- 
tering in the inner crust. 

t = J dra Q (A)n n (r) « Ara (n) (A) (7) 

Here Ar « 500 m is the thickness of the inner crust, 
the nucleon density is n n (r) and (n(r)) as 0.05 fm~ 3 is 
a typical density in the inner crust. Finally (A) is (A) 
appropriately averaged over the composition in the inner 
crust. If the WIMP momentum p is large enough so 
that 2p > q* then a very rough estimate is (A) « 100. 
The mass number A of nuclei present in the inner crust 
tends to increase with density. In addition, A depends 
on the assumed nuclear interaction. For example Thomas 
Fermi calculations of Douchin et al. [l4| or relativistic 
mean field calculations of G. Shen et al. (151 te nd to 
have larger A than semiclassical simulations [16]. We 
will discuss these estimates further in a later paper, see 
also 

With these crude estimates the optical depth of the 
inner crust is r as cto 2.5 x 10 18 fm -2 . Therefore an optical 
depth of order one requires a cross section 

o-q as 4 x 10~ 45 cm 2 . (8) 

Note this is close to Eq. |3J For simplicity we assume 
Mw is (much) less than the mass of the nucleus. In that 



WIMP of momentum p striking a nucleus at rest 
has a maximum momentum transfer of 2p. For WIMPs 
in thermal equilibrium at temperature T one has p as 
(3kTM w ) 1/2 . We require 2p>q*. For M w = 50 GeV 
this is T > 7 x 10 7 K, or 3.5 x 10 8 K for M w = 10 GeV. 
At lower temperatures than these, WIMPs may travel too 
slowly to scatter efficiently in the crust. We conclude that 
it is possible for dark matter to have a significant optical 
depth in the inner crust. However this likely requires 
WIMP cross sections near present experimental bounds. 

WIMPs may be captured, thermalized, and trapped 
using a combination of scattering in the core and in the 
crust. When a WIMP first nears a neutron star, gravity 
accelerates it to a high energy. These energetic WIMPs 
can scatter efficiently in the core because nucleon Pauli 
blocking plays less of a role. The energetic WIMP can 
easily excite nucleons above the Fermi sea so o~wn ~ ctq. 
However, energetic WIMPs can scatter at large q where 
the nuclear form factor F(q) <C 1. Therefore scattering 
in the crust may be inefficient for high energy WIMPs. 
Over time the WIMPs can loose energy by repeatedly 
scattering from nucleons in the core until the WIMPs 
have much lower energy. At that point Pauli blocking 
is important so scattering in the core becomes inefficient 
o~ Wn <C (To- However low energy WIMPs can now scatter 
efficiently in the crust. 

We turn now to transport properties and the contri- 
bution of WIMPs to the shear viscosity of neutron star 
matter. For simplicity we consider just the core of a 
neutron star. We model the core as composed of dense 
neutron rich matter and a low density gas of WIMPs. 
From simple kinetic theory, the shear viscosity of this 
WIMP gas rjw is 

r, w = ^{?,kTM w ) l ' 2 \ w n w . (9) 

Here nw is the number density of WIMPs. 

If the mean free path in Eq. [2] is achieved, let us com- 
pare the shear viscosity in Eq. [9] to the shear viscosity of 
neutron rich matter rj n without any WIMPs. At a den- 
sity of 2no = 0.32 fm -3 and a temperature T « 10 8 K 
(which is typical for a LMXB) Shtcrnin et al. calculate 
rj n = 2 x 10 19 g cm -1 s _1 [17]. This is the sum of elec- 
tron (dominant), muon, and nucleon contributions in a 
non-superfluid star. For Mw = 50 GeV, r)w in Eq. [5] 
will be equal to r\ n if the density of WIMPs is 

n w as 3 x 10 28 cm" 3 . (10) 

This corresponds to a number fraction of WIMPs to nu- 
cleons Xw of 

AV = ^i«10- 10 . (11) 

We conclude that WIMPS will likely dominate the shear 
viscosity of dense neutron rich matter if their number 
fraction Xw is larger than this value. This estimate is 
for the core of the star. We expect qualitatively similar 
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estimates for the inner crust. However this should be ver- 
ified in future work. The density of WIMPs in the core of 
a neutron star may depend on the age of the star, the dis- 
tribution of dark matter in the galaxy, and many WIMP 
properties such as mass, boson/fermion nature, and in- 
teractions with both normal matter and other WIMPs. 
If WIMPS significantly enhance the shear viscosity this 
will increase the damping of collective oscillations such 
as r-modes and may impact the maximum spin rate of 
neutron stars. 

How can neutron stars spin so fast? For example, a 
pulsar in Terzan 5 spins at 716 Hz 18]. This is puz- 
zling because rapidly rotating neutron stars appear to be 
unstable to r-mode oscillations These modes are a 
class of oscillations where the restoring force is the Cori- 
olis force [201 ] . The emission of gravitational waves can 
excite r-modes and cause the amplitude of oscillations 
to grow. Bildsten and Andersson et al. first suggested 
that r-modes could provide a limit on the s pin of neutron 
stars in low mass X-ray binaries (LMXB) [21(. Here an- 
gular momentum gained from accretion is radiated away 
in gravitational waves. This keeps the stars from spin- 
ning faster. Indeed, gravitational waves from r-modes 
may be detectable in large scale interferometers such as 

LIGO m. 

An r-mode is unstable if the rate energy is gained 
through gravitational wave radiation (tapping the stars 
rotational energy) exceeds the rate energy is lost through 
damping due to shear and bulk viscosities. This helps 
motivate many calculations of the shear viscosity of neu- 
tron matter |23l| , neutron rich matter with electrons and 
muons I17l l24| , nuclear pasta with complex non-spherical 
shapes |25l ] , and a crust-core boundary layer [2^, [2?} ■ 
In addition, the bulk viscosity of neutron rich matter 
24 , l28l [29j . hyperon matter [30l - [36| . and quark matter 
37H47| has been calculated. The damping from all of 
these different sources appears to be too small to sta- 
bilize the r-modes in many observed neutron stars [l9| . 
This problem is solved if WIMPs increase the shear vis- 
cosity enough to stabilize the r-modes. 

The viscosity in Eq. [5] depends only weakly on My/ 
and T. However, heavy WIMPs may travel a long time 
between collisions and have a long relaxation time. This 
could imply a frequency dependence for r\y/ if the WIMPs 
can not keep up with an r-mode, whose frequency is of 
order the rotation frequency of the star. This frequency 
dependence of r]y/ may depend more strongly on My/ and 
T and could also lead to a nonzero bulk viscosity. Possi- 



ble WIMP contributions to the bulk viscosity should be 
examined in future work. WIMPs can also contribute to 
the thermal conductivity. This may not be important for 
the core of the star because the thermal conductivity is 
thought to be very high even with out WIMPs. However, 
WIMPs could have an impact on crust cooling after ex- 
tended periods of accretion [48 52] . Here accretion heats 
the crust and then the crust is observed to cool rapidly 
after accretion stops. This rapid cooling suggests a very 
high crust thermal conductivity that is likely consistent 
with (or without) an additional WIMP contribiton. In 
addition WIMPs could transport heat in the crust of a 
strongly magnetized star where it may be otherwise dif- 
ficult to transport heat perpendicular to the magnetic 
field. This could lead to more equal surface tempera- 
tures and less time dependence to the thermal X-ray flux 
as the star rotates. 

We speculate that observations of neutron star rotation 
may provide indirect information on dark matter. For 
example young neutron stars may have captured little 
dark matter and as a result their r-modes may become 
unstable at relatively low frequencies. Therefore one may 
not expect to observe young rapidly rotating stars. In 
contrast old neutron stars may have captured more dark 
matter, stabilizing their r-modes. Indeed a number of old 
neutron stars are observed to spin quickly as millisecond 
pulsars. Alternatively, neutron stars may spin quickly in 
regions of the galaxy with a high dark matter density, 
such as the galactic center. Note that the neutron star 
with the fastest known rotation is near the galactic center 

In summary, neutron stars are attractive places to look 
for dark matter because their high densities allow re- 
peated interactions. Weakly interacting massive parti- 
cles (WIMPs) may scatter efficiently in the core or in 
the crust of a neutron star. In this paper we focused 
on WIMP contributions to transport properties because 
these can be greatly enhanced by long WIMP mean free 
paths. We speculate that WIMPs increase the shear vis- 
cosity of neutron star matter and help stabilize r-mode 
oscillations. As a result, neutron stars can spin rapidly. 
This likely requires WIMP-nucleon cross sections near 
present experimental limits and an appropriate density 
of WIMPs in neutron stars. 
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